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The leading Pollicott-Ruelle resonance is calculated analytically for a general class of two-dimensional area- 
preserving maps. Its wave number dependence determines the normal transport coefficients. In particular, a 
general exact formula for the diffusion coefficient D is derived without any high stochasticity approximation 
and a new effect emerges: The angular evolution can induce fast or slow modes of diffusion even in the high 
stochasticity regime. The behavior of D is examined for three particular cases: (i) the standard map, (ii) a 
sawtooth map, and (iii) a Harper map as an example of a map with nonlinear rotation number. Numerical 
simulations support this formula. 
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Diffusion is a paradigm of deterministic chaos and its study is not new, dating back to Chirikov fl]. Its existence in Hamil- 
tonian systems has been extensively established using a variety of approaches J2, S H IH S 0] • However, it was not clear that 
a satisfactory transport theory could be properly formulated. In order to understand deterministic diffusion, nonequilibrium 
statistical mechanics was suitably combined with dynamical system theory In this modern formulation, the stochastic 

properties of chaotic systems can be determined by the spectral properties of the Perron-Frobenius operator U. One of the most 
important properties is the exponential relaxation to the thermodynamic equilibrium, explained in great detail at the microscopic 
level. The relaxation rates j m , known as Pollicott-Ruelle (PR) resonances iTToL H], are related to the poles z m of the matrix 



elements of the resolvent R(z) = (z — U) 1 as z rn = e 7m . These resonances are located inside the unit circle in the complex 



z plane, whereas the spectrum of U is confined to the unit circle because of unitarity 01211 . Furthermore, the wave number de- 
pendence of the leading PR resonance determines the normal diffusion coefficient for spatially periodic systems These 
results are rigorous only for hyperbolic systems, though they have been confirmed in the high stochasticity approximation for 
some mixed systems such as the kicked rotor (standard map) lfl5ll . the kicked top [16] and the perturbed cat map lfl7ll . The PR 
resonances are essential not only in classical dynamics but also in quantum dynamics. Recently, a microwave billiard experi- 
ment demonstrated a deep connection between quantum properties and classical diffusion through the spectral autocorrelation 
function iflill . 

In this Letter, the leading PR resonance will be calculated analytically for the general class of two-dimensional area-preserving 
maps 

I n+1 =I n +Kf(8 n ), 

dn+i = n + ca(I n+1 ) mod27r, 

defined on the cylinder — ir < 8 < tt, — oo < / < oo. Here f(8) is the impulse function, a(I) = a(I + 2irr) is the rotation 
number, c and r are real parameters, and K is the stochasticity parameter. This map is commonly called the radial twist map lfl9ll 
periodic in action (the nonperiodic case can be considered in the limit r — > oo). Although considerable theoretical development 
in the study of diffusion has been achieved for the linear rotation number (LRN) case ca{I) = I (2j,|3i|4j,|5l|6J], many physically 
realistic systems are best described just by the nonlinear cases. Such maps have been extensively used in various areas of 
physics, especially in celestial mechanics [20], plasma and fluid physics 112 111 , and astrophysics and accelerator devices 
However, the normal transport properties of such maps have not been studied previously |2|]. 

The analysis of the map (Q~|) is best carried out in Fourier space. The Fourier expansion of distribution function at the ?ith time, 
denoted by p n , is given by 

p n (I, 6) = ]T J dqe^ m9+ ^a n (m, q) . (2) 



The moments can be found from the Fourier amplitudes via (I p ) n — (27r) 2 [(i d q ) p a n (q)]q=o, where a n (q) = a n (0, q). The 
discrete time evolution of the probability density p is governed by the Perron-Frobenius operator U defined by p n+ i(I, 9) = 
Up n (I, 9). The matrix representation of U may be considered as the conditional probability density for the transition of the 
initial state (/', 9') to a final state (I, 9) in a one time step, ruled by (Q}. The law of evolution of the Fourier coefficients will be 



2 



given by 



i n (m,q) J dq' A m {r,c,q' - q) J m - m '(-Kq')a n -i{m',q') . 



where ao(m, q) = (2ii) 2 exp[— i(m0 o + qlo)]. The Fourier decompositions of the a(I) and f(6) functions are 

Am(r,c,x) = 7^ S (Zr~ x - x) Gi (r, mc) , 



Gi(r,x) = ^ J d6 exp{-i[xa(r9)-W}}, 
J m {x) = ^[ d6exp{-i[™0-xf(0)]} 



(3) 

(4) 

(5) 
(6) 



If the rotation number a(I) is an odd function, then Gi(r,x) is a real function and G±\i\(r, x) — G\i\{r, ztx). The integral 
function J m {x) assumes the following series expansions J m {x) = S m fi + XmLi c m,n x n , whose coefficients are given by 



1 i n 
2?rn! 



d6f n {6)e 



(7) 



Note that, if f(—6) = —f(9), the coefficients c mj „ are real for all {m, n} and J_ m (x) = J m (—x). 

Let us consider the decomposition method of the resolvent R{z) based on the projection operator techniques utilized in B23I1 . 
The law of evolution J3| can be written as a n (m, q) — U n ao(m, q) where U n is given by the following identity § c dzR(z)z n = 
2iriU n , where the spectrum of U is located inside or on the unit circle C around the origin in complex z-plane. The contour 
of integration is then a circle lying just outside the unit circle. We then introduce a mutually ortoghonal projection operator 
P = \q, 0) {q, 0|, which picks out this relevant state from the resolvent and its complement Q = 1 — P, which projects on 
the irrelevant states. In order to calculate the diffusion coefficient D = lim n -»oo (2w) _1 ((/ — Io) 2 ) n we can decompose the 
projection of the resolvent PR{z) into two parts: PR(z) = PR(z)P + PR(z)Q. The last part can be neglected because 
ao(m y£ 0, q) oc e -jm00 , whose expected value disappear at random initial conditions on [— n, tr). Hence, the relevant law of 
evolution of the Fourier amplitudes, omitting initial angular fluctuations, assumes the following form: 



a n{q) = 77— 

2m 



dz ■ 



c 



where the memory functions 9j (q) obtained for the system (fl]i are given by 

= ^J- m {-Kq)J m {-Kq) Ga(r,mc) 



(8) 



(9a) 



(9b) 



*i>a(?) 



{m} {A}t 



n Jmi (-Kq) G\ t (r, mic) J| G\i (r, rmcjJm^-j 



i=2 



-K 



_1 E^ 

fc=i 



(9c) 



Hereafter, the following convention will be used: Wave numbers denoted by Roman indices can take only non-zero integer 
values, whereas wave numbers denoted by Greek indices can take all integer values, including zero. For each fixed j, the sets 



,} and {A}t = {Ai,...,Aj}, where the superscript denotes the restriction 



of wave numbers are defined by {to} = {mi, . . 

v; , a, o. 

For usual physical situations (assumed here) we have co,i oc J d0f(9) = B2411 . In this case, ^>o(q — > 0) = 1 + 0(q 2 ). 
In the general case, we have >J j(q — > 0) = 0(q 2 ) for j > 1. The integral (O can be solved by the method of residues 
truncating the series at j = N and after taking the limit N — * oo. The trivial resonance z = 1 is related to the equilibrium 
state found for m = m' = q = 0. The nontrivial leading resonance can be evaluated by the well-known Newton-Raphson 
iterative method beginning with z = 1 and converging to = J^Lo ^jio) + ^(l 4 )- ^ n tr,e limit q — > 0, this resonance will 
dominate the integral in the asymptotic limit n — » oo. Thus, the evolution of the relevant Fourier coefficients can be written as 
a n {q) = exp[nj(q)] a (q), where the leading PR resonance is given by 



7 (q)=lnJ2^M+0(q 4 



(10) 



3=0 
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From ([Tol l the diffusion coefficient can be calculated as D = —(1/2) d 2 • Applying this expression to the 

memory functions d9ab-(l9cT>. the general exact diffusion coefficient formula will be given by 

jj oo oo j 

— = 1 + 2^2 o- m , m Re[g (r, mc)} + X! X! X! ^."j Gx t (r, mic) JJ Qx^r, True) l 7m 1 _i-m 4 ( -— X Ai 

9 ' m=l j=2 { m } {A}t i=2 

where D q i = —cq^K 2 is the quasilinear diffusion coefficient and = (c_ TOi i c m ',i)/co,2- The diffusion formula (fTTt 

assumes a more simple form for the LRN case (where / can be replaced by / mod 2ir), yielding Q\ (1, a;) = 5\ tX and 




=2 { m }t 




As a check for this theory, we can first calculate Dlrn explicitly for two cases: (i) the well-known standard map (sm) 
as an example of a mixed system and (ii) a sawtooth map (sw) as a example of a hyperbolic system in a certain parameter 
regime. In case (i), we have f(6) — sin(0), and, hence, J m (x) is the Bessel function of the first kind J m (x), D q i = K 2 /4, 
and cr m ,m' = (±(5 mi ±i)(±i5 m ' i ±i). The resultant expression for D sm is very similar to ( fl~2b . The first terms of the expansion 
coincide with the Rechester, Rosenbluth and White results [2]: D sm /D q i = 1 — 2J2(K) + 2Jf (AT) + In case (ii), we 

have f(8) = 8; hence, J m {x) — ~fe^z > -Dgi — K 2 n 2 /6, and Cm.m' = -J^ ''"mm' — ■ ^h e sawtooth ma P is hyperbolic 
when |A' + 2| > 2 N25I1 . Finally, we also consider a third case where f(x) = a(x) — sin(a;), known as the Harper map 
(Hm), as an example of a map with a nonlinear rotation number. The resultant expression for Dnm is very similar to (fTTT i. 
where <5a(1, a;) = J\{x) and other terms follow case (i). The analytical results on the three cases are compared with numerical 
calculations of D/D q i in Figs. l(a)-l(c). Despite the accelerator modes, whose kinetics properties are anomalous ll26ll . all 
theoretical results are in excellent agreement with the numerical simulations. 

A question of interest that arises here is the oscillatory character of the diffusion coefficient for maps with a periodic rotation 
number (including the LRN case), in contrast to the fast asymptotic behavior exhibited by maps with a nonperiodic rotation 
number (see, for example, J3l). The nonperiodic case can be considered by applying the limit r — > oo. For such a case, we 
have Ar _1 — > s, r -1 ~J2 X — > J ds, and r 0\(r, x) — > Q(s,x) is the s-Fourier transform of e~ ixa ^ I \ For cases where Q(s,x) 
is well defined, the limit r — > oo produces high oscillatory integrals resulting in D — > D q i without any oscillation. In the case 
of a standard map, Chirikov QjJ] conjectured that the oscillatory aspect of the diffusion curve was an effect of the "islands of 
stability", but a satisfactory explanation for the oscillations has not been given yet 12711 . 

Returning to equation (fTTT > we can note that, in the limit of high stochasticity parameter K, the diffusion coefficient does not 
necessarily converge to the quasilinear value in the nonlinear rotation number cases. The standard argument in this respect is the 
so-called random phase approximation QjJ, Ll9|] - The intuitive idea is that, for large values of K, the phases 9 n (I, 8) oscillate so 
fast that they become uncorrected from 8. In order to verify this effect we can take the limit K — * oo of ( fTTT i by setting \ — 
for all i to avoid terms of order 0{K~ X ' 2 ). Once \Go{r, mc) \ < 1 for m ^ 0, the asymptotic difusion becomes a geometric sum 
whose result is 

Jim — = 1 + > <J m , m - tt-, r- (13) 

K^oo D q i *-f l-g (r,mc) 

The rate ( fT3l diverges at c = 0, creating a kind of accelerator mode. Indeed, a direct calculation through Eq.(fl]i shows that 
D/D q i diverges as n in this case for all K ^ 0. In Fig. 1(d) we consider the double sine map for K — 10 5 . As one can 
see, even in the high stochasticity regime, where the random phase approximation is expected to hold, the rate D/D q i oscillates 
between the zeros of Jo(c). Its maximum and minimum values are ruled by zeros of J\ (c). This strong angular memory effect 
is a remarkable result. 

Another important question concerns the higher-order transport coefficients that play a central role in the large deviations 
theory. These coefficients can be obtained through the following dispersion relation: 

^■J5i i5 W 2 *-W fl!S( " u ' <14) 

where I > 1 and ( } c denotes cumulant moments [28]. The diffusion coefficient is defined by D = T>2- The higher-order coeffi- 
cients D21 can be calculated by introducing successive corrections 0(q 21 ) in (TTOb . If the evolution process were asymptotically 
truly diffusive, then the angle-averaged density wold have a Gaussian contour after a sufficiently long time. A first indication of 
the deviation of a density function from a Gaussian packet is given by the fourth-order Burnett coefficient B = V4: If B = 0, 
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then the kurtosis k(x) — (x 4 ) / ^x 2 ) for x — I n — Iq is equal to 3 in the limit n — > oo, a result valid for a Gaussian density for 
all times. These aspects will be treated elsewhere ll29ll . 

The author thanks Professor A. Saa, Professor W.F. Wreszinski, and Professor M.S. Hussein for helpful discussions. This 
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FIG. 1: Theoretical diffusion coefficient rate D/D q i (solid lines) compared with numerical simulations calculated for n = 100. In the cases 
(a), (b), and (c) we truncate the diffusion formulas i ll It and d!2t at j — 2. A better agreement for small values of K requires the calculation 
of further memory functions, (a) Standard map. The accelerator modes give rise to spikes in the figure, (b) Sawtooth map and (c) Harper map 
for c = 5.5 (with presence of accelerator modes), (d) Harper map as a function of c for K = 10°. The angular evolution induces fast and slow 
modes of diffusion even in the high stochasticity regime. This strong angular memory effect decay as 2Jo(c)/[l — Jo(c)]. 



